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^ : ABSTRACT 
(N ■ 
(N ■ 

\ The escape fraction, /esc, of ionizing photons from early galaxies is a crucial 

Q \ parameter for determining whether the observed galaxies at 2; > 6 are able to 

^ \ reionize the high-redshift intergalactic medium. Previous attempts to measure 

\ /esc have found a wide range of values, varying from less than 0.01 to nearly 
1. Rather than finding a single value of fesc, we clarify through modeling how 

O I internal properties of galaxies affect /esc through the density and distribution 

I of neutral hydrogen within the galaxy, along with the rate of ionizing photons 

I production. We find that the escape fraction depends sensitively on the covering 

^ I factor of clumps, along with the density of the clumped and interclump medium. 

>• ' One must therefore be cautious when dealing with an inhomogeneous medium. 

^ I Fewer, high-density clumps lead to a greater escape fraction than more numerous 

^ I low-density clumps. When more ionizing photons are produced in a starburst, /esc 

^ ' increases, as photons escape more readily from the gas layers. Large variations in 

O ' the predicted escape fraction, caused by differences in the hydrogen distribution, 

^ ■ may explain the large observed differences in /esc among galaxies. Values of /esc 

' must also be consistent with the reionization history. High-mass galaxies alone 

■ are unable to reionize the universe, because /esc > 1 would be required. Small 

;h ' galaxies are needed to achieve reionization, with greater mean escape fraction in 

■ - - ' the past. 
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INTRODUCTION 



Observations of the cosmic microwave 



jackground optica . 



son Microwave Anisotropy Probe fWMAP) (IKogut et al 



20071 : ISpergel et all 120071 : iDunklev et al.ll2009l : iKomatsu et al. 



2003 



depth made with t he Wilkin- 



Sperge 



2009 



et al.ll2003l : IPage et al 



20101) suggest that the 



universe was reionized sometime between 6 < z < 12. Because massive stars are efficient 
producers of ultraviolet photons, they are the most likely candidates for the majority of 
reionization. However, in order for early star-forming galaxies to reionize the universe, their 
ionizing radiation must be able to escape from the halos, in which neutral hydrogen (H i) is 
the dominant source of Lyman continuum (LyC) opacity. The escape fraction, /esc, of ioniz- 
ing photons is a key parameter for starburst galaxi es at 2: > 6, which are believed to produce 



the bulk of the photo ns that reionize the universe (IRobertson et al.ll2010l : iTrenti et al.l 12010 
Bouwens et al.ll2010l ). 



The predicted values of escape fraction span a large range from 0.01 < /esc < 1, derived 
from a variety of theoretical and observational studies of varying complexity. Various proper- 
ties of the host galaxy, its stars, or its environment are thought to affect the number of ioniz- 
i ng ph otons that escape into the intergalactic medium (IGM). For ex ample. iRicotti fc Shull 
(120001 ) studied /esc in spherical halos using a Stromgren approach. IWood fc Loebl J200oh 
assumed an isothermal, exponential disk galaxy and followed an ioniza tion front through the 
gala xy using three-d i mensi onal Monte Carlo radiative transfer. Both IWood &: Loebl (120001 ) 
and iRicotti fc Shulll (I2OOOI ) state that /esc varies greatly, from < 0.01 to 1, depending on 
galaxy mass, with larger galaxies giving small er values of f e sc A similar dependence with 



galaxy mass is also seen by the simulations of lYajima et al.l (120101 ). because larger galaxies 
tend to have star formation buried within dense hy drogen clouds, while smaller galaxies often 
had clearer paths for escaping ionizing radiation. iGnedin et al.l (120081 ) . on the other hand, 
ran a high resolution N-body simulation w i th adaptive-mesh ref inem ent in a cosm o logica l 
context. Contrary to IRicotti fc Shulll J200oh . lWood fc Loebl J200oh and lVaiima et all toid ). 
they state that lower-mass galaxies have significantly smaller /esc, as the result of a declining 
star formatio n rate. In add i tion, a bove a critical halo mass, /esc does not change by much. 
The model of iGnedin et al.l (120081 ) allowed for the star formation rate to increase with the 
mass of the galaxy at a higher rate than a linear proportionality would allow. The larger 
galaxies also tended to have star formation occurring in the outskirts of the galaxy, which 
made it easier for ionizing photons to escape. Their model inclu ded a distribution of gas 
within the galaxy, which created free sight-lines out of the galaxy. IWise Sz CenI (120091 ) used 
adaptive mesh hydrodynamical simulations on dwarf galaxies. Even though their s imula - 
tions covered a different mass range than the larger galaxies studied by IGnedin et al.l (120081 ). 
they found much higher value of /esc than would be expected from extrapolating results 
from IGnedin et al.l (120081 ) to lower masses. IWise fc CenI (120091 ) attribute this difference to 
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the irregular morphology of their dwarf galaxies with a turbulent and clumpy interstellar 
medium (ISM), allowing for large values of /esc- 

Others have als o looked at how the shape and morphology of the galaxy can affect /esc- 
Dove fc Shulll (119941 ). using a Stromgren model, studied how /esc varies with various H i disk 
density distributions. In addition, many authors have foun d that supe r bubbl es and shells 
can trap radiation until blowou t, seen in ana l ytical models of iDove et al.l fl2000l) as well as in 



hydrq dynamical simulations of iFujita et al.l ( l2003l ). The analytical model by IClark &: Oey 
( I2OO2I ) showed that high star formation rates can raise the porosity of the ISM and thereby 
increase /esc- In addition to bubbles and structure caused from supernovae, galaxies can 
have a clur npy ISM whose i nhomogeneities affect fmc- For ex a mple, dense clump s coul d 



reduce f^^r- (iDoye et a. 



2000). On the othe r hand. iBoissel f ll990h . iHobson fc Scheued ( 119931 ). 



Witt fc GordonI ( Il996h . and lWood fc Loebl feoooh all found that clumps in a randomly dis- 
tributed medium cause /esc to rise, while ICiardi et al.l ( 12002| ) found that the effects of clumps 
depend on the ionization rate. 

A host of other galaxy parameters have been tested analytically and with simulations. 
Increasing the baryon mass fraction lowe rs /esc for smaller halos, but increases it at masses 
greater than 10^ Mq (IWise fc CenI 120091 ). Star formation history changes the amount of 
ionizing photons and neutral hydrogen, causing /esc to vary from .12 to 0.20 for co eval star 
formation and from 0.04 to 0.10 for a time-distrib uted starburst (iDove et al.ll2000l). O ther 
galactic quantities, such as spin (IWise fc Cenll2009l ) or dust content ( iGnedin et al.ll2008l ). do 
not seem to affect the escape fraction. 

Observations have also been used to constrain /esc, especially at 2; < 3. Searches for 
escaping Lyman co ntinuum radiation at redshifts z 5, 1 — 2 have found escape fractions of at 



most a few percent (Bland-Hawthorn fc Malonev 



ll999[ Deharveng. et al. 



2002 



Bridge et al. 



i 



2010 



Cowie et al. 



2009 



Tumlinson et a 



Leitherer et al 



19951: 



Malkan. et al. 



2001 



Grimes et al. 



2003[ISiana. et al.lboOTl) . 



variations in the escape fraction, and iHoopes et al.l ( 2007 ) and iBergvall et al.l (|2006[ ) saw a 
relatively high esca pe fraction of 10%. iFergusonI (l200l[ ) observed /esc ~ 0.2 at 2; ~ 1. 
Hanish et al.l (I2OIOI) do not see a difference in /esc between starbursts and normal galaxies. 



200 



20091: 



Heckman et al. 



Hurwitz et al 



2001 



l|l997|l saw large 



Siana et al. ~( l2010h also found low escape fractions at z ~ 1.3 and showed that no more than 
8% of galaxies at this redshift can have /esc,rei > 0.5. Note that /esc.rei, which the authors use 
to compare their results to other surveys, is defin ed as the ratio of escaping L y C ph otons 
to escaping 1500 A photons. In our own Galaxy, iBland-Hawthorn fc Maloneyl (1999) and 



Putman et al. 



(I2OO3I) found an escape fraction of only a few percent. Observations using 7 
ray bursts (IChen et al.ll2007l ) show /esc ~ 0. 02 a.t z ^ 2. At higher redshift (z 3), fesc seems 



to vary drastically from galaxy to galaxy ( IShapley et al.l l2006l : llwatal 120091 : IVanzella et al. 
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2010bl ). with a few galaxies having very larg e escape fractions. Some studies have found 



low values of the escape 



2002: Heckman et al. 2001 



significant LyC leakage ( ISteidel et al. 



Taction at z 



Inoue et a. 



2005 



2002 



Wvithe et a. 



Shapley et al. 



boioh 



Fernandez- Soto et al.l l2003l : iGiallongo. et al. 



.. ■■ while others have found 
20061 ). This large variation from 



galaxy to galaxy suggests a dependence on viewing a ngle and could indicate the patchiness 



and structure of neutral hydrogen within the g alaxy (jBergvall et al.ll2006l : iDeharveng. et al. 



200ll : iGrimes et al.ll2007l : iHeckman et al.ll200ll ). From these observations, one infers that the 



fundament al properties o 



the galaxies charige with time or that /psr increases with increas- 



ing red shift flBridge et al.ll2010l : iCowie et al.ll2009l : llnoue et al.ll2006l : Ilwatall2009l : iSiana. et al. 
20071 ). iBouwens et al.l (120 lOf ) looked at the blue color of high redshift (z ~ 7) galaxies and 
argued that the nebular component must be reduced. This would suggest a much larger 
escape fraction in the past. 

The minimum mass of galaxy formation can also put limitations on /esc- Observations of 
Lya absorption toward high-redshift quasars, combined with the UV luminosity function of 
galaxies, can limit /esc from a redshift of 5.5 to 6, with /esc ~ 0.20-0.45 if the halos producing 
these photons are larger than 10^° Mq. This can decrease to /epc ~ 0.05-0.1 if halos down to 



10^ Mq are included as sources of escaping ionizing photons (jSrbinovsky fc Wyithdl2008l ). 



It is clear that many factors can affect /esc, and the problem is quite complicated. Cos- 
mological simulations that predict the escape fraction provide a more accurate estimate for 
/esc- However, many parameters of the galaxy change at once, and it becomes difficult to 
understand how a single parameter can affect the escape fraction. In addition, trends may 
be difficult to understand because of the manner in which some physical processes are in- 
cluded or neglected. Analytic models can show clearer trends, even though they may be 
over-simplified and miss important physics. Therefore, rather than predicting a quantitative 
value for /esc, we seek to understand how properties of galaxies and their internal structure 
affect the escape fraction. Because our model is simplified, the values of /esc are not exact, 
but rather illustrate trends caused by various galactic properties. In section [2], we explain 
our method of tracing photons that escape the galaxy. In section [3], we explain our results 
and compare our results to previous literature in section |H In section [5|, we consider con- 
straints from reionization an d we conclude in secti on El Throughout, we use the cosmological 
parameters from WMAP-7 (IKomatsu et al.ll2010l ). 
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METHODOLOGY 



2.1. Properties of the Galaxy 



We use an exponential hyperbolic secant profile (ISpitzerlll942l ) to describe the density 
of an isothermal disk in a halo of mass Mhaio: 



nniZ) = no exp[— r/r/j] sech^ 



( lSpitzerlll942[ ) where rio is the number density of hydrogen at the center of the galaxy, Z is 
the height above the galaxy mid-plane, and r/, is the scale radius: 



(2) 



( IMo et al.lll998l ). The parameter jd is the fraction of the halo's angular momentum in the 
disk, A is the spin parameter, r rid is the fracti o n of t he halo in the disk (m^ = Qb/^m), and 
Vyir is the virial radius. As in IWood fc Loebl (120001 ). we assume jd/^d = 1 and A = 0.05. 
The virial radius is 

i?xi;)2oo' ' — ' ' 



(0.76kpc) 



V lO^Mo/i 



10 



(3) 



(INavarro et al.lll997l ) where Ac = ISvr^ + 82d — 39d^ and d = Qzf — 1- ^zf is the local value 
of flm at the redshift of galaxy formation, zj. The dependence of the virial radius on Zf will 
affect the density of the disk, with smaller disks of higher density forming earlier. The disk 
scale height, zq, is given by 



\27rGpo 



1/2 



M, 



halo 



1/2 



(4) 



V27rpo^^. 

where (f ^) is the mean square of the velocity and po is the central density. In a real galaxy, 
there will be non-thermal motions of the gas. However, to simplify the calculations, we 
assume that the gas is virialized, feels the gravity of the disk, and therefore follows the 
relation (u^) = GMhaio/fvir- However, radiative cooling can cause the disk to be thinner 
than this. The central density is solved for in a self-consistent way after the halo mass and 
the redshift of formation are specified. We use ISr/^ and 2zq as the limits of the radius and 
height of the disk, respectively. The structure of real galaxies is more complicated. The 
addition of stars in the halo will change the gas distribution through feedback, heating, and 
gravitational effects. For purposes of simplicity, we ignore these effects. 

The mass of the disk (stars and gas) is taken as Mjisk = m^Mhaio; where rrid is the 
fraction of matter that is incorporated into the disk. The upper limit of rrid is ^Ib/^m- The 
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mass of the stars within the disk is M* = Mdisk/*, where /* is the star formation efficiency, 
which describes the fraction of baryons that form into sta;^^ The remainder of the mass of 
the disk is in gas, distributed according to equation [T] with gas temperature of 10"^ K. 

The number of ionizing photons is related to /*, considering either Population III (metal- 
free) or Population II (metal-poor, Z = 0.02Zq) stars. The total number of ionizing photons 
per second from the entire stellar population, Qpop is 

C QHijn)f{rn)dm 



Q 



pop 



mf{m)dm 



X 



(5) 



where m is the mass of the star, and mi and m2 are the upper and lower mass limits of the 
mass spectrum, given by ffm). For a less massive distribution of stars, we use the Salpeter 
initial mass spectrum (ISalpeterlll955[ ): 



/(m) oc m 



-2.35 



(6) 



with nil = 0.4 Mfr^ and mo = 15 Mq. The Larson initial mass spectrum illustrates a case 
with heavier stars (ILarsonlllQQSl ): 

-1.35 

(7) 



m 

1 + — 

mr 



f{m) (X m ^ 

500 Mq, and m^ = 250 Mq for Population III stars and mi 



1 M 



with mi = 1 Mq, m2 

7712 = 150 Mq, and mc = 50 Mq for Population II stars. We define Q^j as the number 
of ionizing photons emitted per second per star, averaged over the star's lifetime. For 
Population III stars of mass parameter x = logiQ^m / Mq) , this is 



logio \Qh/s ^] 
and for Population II stars, 
logio \Qh/^~^] = 




4.90x - 0.83x2 g _ ^ 



1 



8.55x 



5 - 9 M© , 
otherwise , 



27.80 




30.68X - 14.80x2 + 2.50x3 > 5Mq 

otherwise 



(9) 



as given in Table 6 of ISchaererl (120021 ). 



^The star formation efficiency is defined by tlie fraction of baryons in stars at any given time. Therefore, 
the escape fraction calculated is the escape fraction at that point in the galaxy's lifetime. The total escape 
fraction of the galaxy will depend on the lifetime of the burst and the star formation history, which can be 
obtained by integrating over the duration of the burst, taking into account the propagation of the ionizing 
front within the galaxy. 
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2.2. Calculating the Escape Fraction 

We place the stars at the center of the galaxy. An ionized H II region develops around 
the stars, where the number of ionizing photons emitted per second by the stellar population, 
Qpop, is balanced by recombinations, such that 

Qpop = ^TrryHaBiT). (10) 

Here, is the case-B recombination rate coefficient of hydrogen and T is the temperature 
of the gas (we assume T = W^K). The radius of this H ii region, called the Stromgren 
radius, is 



3Qpop ^ 



This radius is simple to evaluate in the case of a uniform medium, but if we are concerned 
with clumps and a disk with a density profile, the density will be changing with location. 
Although we are calculating /esc at a moment in time, in reality the H ii region is not static, 
and the ionization front will propagate at a fiux-limited speed. 

We assume that all stars are placed at the center of the galaxy. In reality, star formation 
will be distributed throughout the galaxy. If stars are closer to the edge of the galaxy, their 
photons will have less hydrogen to traverse, and hence will escape more easily. Therefore, 
the results we present will be lower limits of the escape fraction. We integrate along the 
path length t h at a p hoton takes in order to escape the galaxy, following the formalism in 



Dove fc Shulll (119941 ). We can then calculate the escape fraction of ionizing photons along 
each ray emanating from the center of the galaxy by equating the number of ionizing photons 
to the number of hydrogen atoms across its path. If there are more photons than hydrogen 
atoms, the ray can break out of the disk; otherwise, no photons escape and the escape fraction 
is zero. The escape fraction along a path, rj, thus depends on the amount of hydrogen the 
ray transverses, which depends on its angle 6, measured from the axis perpendicular to the 
disk: 

V{0) = 1 - / nUzydr . (12) 

^pop J 

Photons are more likely to escape out of the top and bottom of the disk, rather than the 
sides, because there is less path length to traverse. This creates a critical angle, beyond 
which photons no longer will escape the galaxy. The total escape fraction, /esc, is then found 
by integrating over all angles 6 and the solid angle Q: 

/esc(gpop) = // "^dOdn (13) 

^7]{e) sm{e) de . (14) 
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We take into account the whole disk (top and bottom) so that an /esc of 1 means that all 
photons produced are escaping into the IGM. 



2.3. Adding Clumps 

A medium with clumps can be described with the density contrast C = Uc/nic between 
the clumps (density Uc), and the interclump medium (density Uic). The percentage of volume 
taken up by the clumps is described by the volume filling factor fy- We randomly distribute 
clumps throughout the galaxy. We define nmcan as the density the medium would have if it 
was not clumpy, given by equation [H The density at each point is given by 

fv + {l- fv)/C ^''^ 

if the point is in a clump and 

^mean /-i /^\ 

= MC - 1) + 1 <i«' 



if the point is not in a clump, similar to (jWood fc Loebl 120001 ). In this way, the galaxy 



retains the same interstellar gas mass, independent of fv and C. As C increases, the density 
of the clumps increases as the density of the non-clumped medium falls. Similarly, if fv is 
larger, more of the medium is contained in less dense clumps. We trace photons on their 
path through the galaxy and track whether or not they encounter a clump. At each step on 
the path out of the galaxy, a random number is generated. A clump exists if this random 
number is less than the volume filling factor. As the filling factor increases, clumps can 
merge, forming larger, arbitrarily shaped clumps. Counting the number of photons exiting 
the galaxy then leads to /esc- The covering factor is computed by counting how many clumps 
intersect a ray as it travels out of the galaxy. 



3. RESULTS 

3.1. Properties of the Clumps 

In the first calculations, we placed Population III stars with a Larson mass spectrum 
and a star formation efficiency /* = 0.5 in a halo of Mhaio = IO^Mq, with a redshift of 
formation of Zf = 10. The clumps have diameter 10^^ cm, unless otherwise stated. The 
top panel of Figure [1] shows /esc as a function of fv for various values of C. The case with 
no clumps is equivalent to C = L As clumps are introduced, /esc quickly falls, but rises 
again as fv rises. This is because the clumps become less dense (since more of the medium 
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is in clumps and the mass of the galaxy must be kept constant). In addition, /esc drops as 
C increases, showing that denser clumps with a less dense interclump medium stops more 
ionizing radiation than a more evenly distributed medium. The clumps are small enough 
that essentially every ray traversing the galaxy encounters one of these very dense clumps 
and is diminished. 

In the bottom panel of Figured! the same population of stars is shown for various values 
of fv as a function of log(C). As C increases, /esc becomes low for small values of fy. Again, 
this is because of a few very dense clumps that stop essentially all radiation. As fy increases, 
more of the medium is in clumps, and therefore the density of the clumps decreases. The 
combined effect is an increase of /esc- The solid black line shows the case with no clumps, or 
when fy = 0. For fy = 0, /esc equals the case with no clumps (C = 1), as it should. Above 
C ~ 10 — 100, increasing C no longer affects /esc- 

So far, we have only been exploring the results of small clumps (10^^ cm, or ~ 0.3 pc) in 
diameter. What would happen if we were to increase the size of these clumps? In this case, 
a ray would traverse fewer clumps as it travels out of the galaxy (the covering factor will 
fall), but any given clump would be larger. As shown in Figure [21 fesc rises as the clumps 
increase in size. For very low values of fy, only a few clumps exist and not every ray comes 
in contact with a clump, increasing the escape fraction above the case with no clumps. 

To illustrate this further. Figure E] shows /esc against fy for a large clump size (10^^ cm 
in diameter). In the top panel, C is varied for = 0.01 and = 0.17. The left-most 
vertical line represents the value of fy needed for a photon traversing the longest path length 
to pass through an average of one clump for both = 0.01 and rrid = 0.17. The right-most 
vertical solid line represents the value of fy for a photon traversing the shortest path length 
to pass through an average of one clump for = 0.17, while the dashed line is the shortest 
path for rrid = 0.01. To the right of these lines, all path lengths intersect a clump. In other 
words, here the covering factor is greater than one. To the left of these lines, there are 
clump-free path lengths out of the galaxy. For = 0.17, we see that if fy is low enough 
some rays will pass through fewer than one clump, on average, /esc is much greater than 
the no-clump case. For very low values of fy, there are so few clumps that the interclump 
medium approaches the case with no clumps. Therefore, the plot of /esc is peaked in the 
region where there are some paths that do not intersect a clump. These results are averaged 
over ten runs, while Figure H] shows the distribution for each run. 

The bottom panel of Figure |3] shows /esc for various values of m^. As increases, the 
disk is less massive, and the escape fraction for a non-clumped galaxy rises. For = 0.01, 
the escape fraction of a non-clumped galaxy is ~ 1. In this case, once clumps are added, the 
escape fraction decreases because regions dense enough to stop ionizing radiation are finally 
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introduced. When the covering factor is greater than one, the clumps grow less dense, 
causing the escape fraction to rise again. For all other cases, /esc peaks when the covering 
factor is less than one, and falls below the escape fraction for the non-clumped case when 
the covering factor is greater than one. 

Star formation is more likely to take place in clumps. Star formation will also have 
an effect on gas clumping produced by stellar winds and supernova shells. Because of this, 
clumps are likely to be distributed around locations of star formation. To analyze this effect, 
we have defined a region 100 pc from the center of star formation in the galaxy. Inside 
this region, the volume filling factor is /v,near, and outside of this region, the volume filling 
factor is /y.far- For a case where clumps are only located near to the star formation center, 
/v,far = 0. If /v,far < /v.near, there are fewer clumps in the outer portions of the galaxy than 
near the star formation center. 

Results are shown in Figure [51 The black solid line shows the case with no clumps, where 
/y = throughout, and the red triple-dot dashed line represents the case where fv = 0.3 
throughout. For the case where /v,near = 0.3 and /v,far = 0, clumps are only near stars. This 
results in a higher escape fraction than if the entire galaxy had fy = 0.3. When there are 
some clumps far from star formation, in the case where /v,ncar = 0.3 and /v,far = 0.1, the 
escape fraction falls below the case where if the entire galaxy had fv = 0.3. This is a result 
of clumps becoming more dense as fv increases. 

We have chosen a 100 pc radius of the region near to the star formation. If this region is 
much larger, ~ 1 kpc, the escape fraction will be equal to the case fv = /v.near, where fv is 
the volume filling factor if the galaxy has the same value throughout. If the region is much 
smaller, ~ 10 pc, the escape fraction will be equal to the case where fv = /v,far- Overall, 
the distribution of clumps does not change the escape fraction significantly. 

3.2. Properties of Stars and the Galcixy 

In Figure [6], we analyze how the stellar population affects the escape fraction. In the 
top panel, /* is held constant as fv increases. In the bottom panel, fv is held constant as 
increases. Both plots show metal-free (Population III) stars and metal-poor (Population 
II) stars, as well as stars with a heavy Larson initial mass spectrum and a light Salpeter 
initial mass spectrum. In both cases, /esc is proportional to the number of ionizing photons 
that are emitted by the stars, with heavier stars or stars with fewer metals more likely to 
produce photons that can escape the nebula. This is because when more ionizing photons 
are produced, the critical angle where photons can break free from the halo increases, and 
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hence more photons escape. When the galaxy forms more stars (higher /*), there is the 
added effect of less hydrogen remaining in the galaxy to absorb ionizing photons. Therefore, 
/esc increases greatly as /* increases. 

The star formation redshift, Zf, is varied in Figure [71 As 2;/ increases, the galaxy is 
smaller and more concentrated. Therefore, it is easiest for photons to escape from less 
dense disks at low redshifts. At redshifts where we expect reionization to take place, it is 
harder for photons to escape the galaxy. This problem may be remedied by the high-redshift 
expectation of more massive or metal-free stars with higher values of f^,. 



4. COMPARISON TO PREVIOUS LITERATURE 



As noted in the introduction, there have been many previous studies that calculated the 
number of ionizing photons emitted from high redshift halos, resulting in a wide range of 
values for the escape fraction. Various factors are proposed that affect the num ber of ionizing 
pho tons that escape into t he IGM, in particular the effects of a clumpy ISM. Isoisse ( 1990 ) 
andlWitt fc GordonI fll996[ ) found that clumps increase transmission, and iHobson fc Scheuer 
( 119931 ) found that a three-phase medium (clumps grouped together, rather than randomly 
dis tributed) furth e r incr eases transmission. Very dense clumps (with C = 10^) were studied 
by lWood fc Loebl ( 120001 ). who found that clumps increase /esc over the case with no clumps. 
For very small values of fv, their /esc was very high, because most of the density is in a few 
very dense clumps, and most lines of sight do not encounter a clump. Their clump size is 
13.2 pc, which is similar to our largest clump size, 5 x 10^^ cm. Their results are consistent 
with our findings for clumps with large radius, low fy, and high C, where most rays do not 
encounter a clump. 



Ciardi et al.l ( 120021 ) included the effect of clumps using a fractal distribution of the ISM 
with C = 4 — 8. They noted that this distribution of clumps increases /esc in cases with lower 
ionization rate because there are clearer sight lines. They found that /esc is more sensitive 
to the gas distribution than to the stellar distribution. 



Dove et al.l ( I2OOOI ) reported that /esc decreases as clumps are added. This results from 



the fact that adding clumps does not change the density of the interclump medium. In their 
model, as clumps are added, the mass of hydrogen in the galaxy increases. On the other 
hand, our current method decreases the density of the interclump medium as clumps are 
added or become denser to keep the overall mass of the galaxy constant. Photons are more 
likely to escape along paths with lower density, as in irregular galaxies and along certain lines 
of sight (iGnedin et al.ll2008t IWise fc Cenll2009l ). Shells, such as those created by supernova 



-VI- 



remnants (SNRs), can trap ionizing photons until the bubble blows out of the disk, allowing; 
photons a clear pa th to escape and causing /esc to rise (IDove fc Shulllll994t iDove et al.ll2000 



Fujita et al.ll2003l ). These SNRs or s uperbubbles create porosity in the ISM, and above a 
critical star formation rate, /esc rises (jClark fc Oeyll2002l ). This is similar to what is seen in 
our results. As with a dense clump, a shell will essentially stop all radiation, while a clear 
path, similar to the case with a low /y, allows many free paths along which radiation can 
escape. Our model extends the previous work by enlarging the parameter space. We see 
how low values of /y, changing the clump size, and location of the clumps affect the results. 
This heavy dependence on how the clumping can affect the escape fraction can explain why 
such a variation is seen in observations in the escape fraction from galaxy to galaxy. 

Previous w o rk di ffers as to whether or not /esc increases or decreases with redshift. 
Ricotti fc Shud fl2000[ ) state that /esc decreases with increasing redshift for a fixed halo 
mass, but is consistent with higher escape fraction from dwarf galaxies. This assumes that 
the star formation efficiency is proportional to the baryonic content in a galaxy. However, 
ot her studies seem to indicate that fes c increases with redshift. High resolution simulations 
of iRazoumov fc Sommer-Larseru (120061 ) state that /esc increases with redshift from z = 2.39 
(where /esc = 0.01-0.02) to z = 3.8 (where fesc = 0.06-0.1). This is a result of higher gas 
clumping and lower star formation rates at lower redshifts, causing the e scape fraction to fall. 
At higher redshifts, the simulations of iRazoumov Sz Sommer-Larseru (120101 ) see this trend 



continue, with an f^sc ~ 0.8 at 2; 10.4 that dec 
seen qbseryationally from < z < 7 (jSridge et "aL 



ines with time. This trend has also been 



Cowie et al. 


2009; 


Inoue et a 


2006 


owever. 


Vanzella et al. 


( 


2010al 


point 



out that observational measurements of the escape fraction can be co ntaminated by lower 
redshift interlopers.) Other simulations have given different results. iGnedin et al.l (120081 ) 
say that /esc changes little from 3 < 2 < 9, always being about 0.01-0.03. This difference 
could po ssibly be attributed t o how the models deal with star formation efficiencies within 
galaxies. IWood &: Loebl (I2OOOI ) state that since disk density increases with redshift, /esc will 
fall as the formation redshift increases, ranging from fesc = 0.01-1. We found that fesc 
decreases with increasing redshift of formation, since disks are more dense. However, this 
assumes that the types of stars and /* remain constant. If /* is larger at higher redshifts 
and if stars are more massive and have a lower metallicity (likely), the number of ionizing 
photons should increase, which could cause /esc to increase, despite a denser disk. It is also 
possible that the disk morphology does not exist at high redshifts. Therefore, in order to 
understand the evolution of /esc with redshift, one must understand the evolution of other 
properties, namely, /*, Qpop, and the density distribution of gas within a galaxy. 
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CONSTRAINTS FROM REIONIZATION 



If galaxies are responsible for keeping the universe reionized, there must be a minimum 
number of photons that can escape these galaxies to be consistent with reionization. The 
star formation rate (p) that corresponds to a star formation efficiency /* is given by: 



p{z) = 0.536 Mq yr~^ Mpc' 



0.1 0.02 / I 0.14 



2\ 1/2 



1+Z 

10 



3/2 



(17) 



f|Fernandez fc Komatsull2006l ). assuming a Press-Schechter mass function (jPress fc Schechter 
19741 ). with 

Here, dc is the overdensity, cr(M) is the present-day rms amplitude of mass fluctuations, and 
Mmin is the minimum mass of halos that create stars. 

Similarly, the /esc needed to reionize the universe can be related to the critical star formation 
rate, pcrit, needed to keep the universe ionized: 



Pcrit(2) = (0.012 Mo yr-^ Mpc" 



l + z 



1 3 



Ch/5 

f esc/0.5 ^ 



0.004 

QhyC 



-0.845 



(19) 



( iMadau et al.lll999t IShull fc Trentill2010l ). which results from the number of photoionizations 
needed to balance recombinations to keep the universe ionized at 2 = 7. Here, Ch is the 
clumping of the IGM (which we scale to a typical value of Ch = 5), T4 is the temperature of 
the IGM in units of 10^ K, and QhyC is the conversion factor from p{z) to the total number 
of Lyman continuum photons produced per Mq of star formation. 



QhyC = 



A^Lyc/10 
Pcrit'^r 



63 



(20) 



where A^LyC is the number of Lyman continuum (LyC) photons produced by a star and tree 
is the hydrogen recombination timescale. We assume that Q^yc = 0.004, which is reasonable 
for a low-metallicity stellar population. 

By requiring the star formation rate to be at least as large as the critical value, we can 
solve for the value of /esc needed to reionize the universe. Results are shown in the top panel 
of Figure [8] plotted at two redshifts, z = 7 and z = 10. We also consider stars forming in 
galaxies down to a minimum mass 



(10^ Mg 



+ 
10 



H -1.5 



(21) 
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( iBarkana fc Loebl l200ll ). or if smaller halos are suppressed and only those above lO^M© 
are forming stars. If the required /esc exceeds 1 (shown by the dashed line), the given 
population cannot reionize the universe. As redshift decreases, it becomes much easier to 
keep the universe reionized, and a smaller /esc is needed, as expected. 

At higher redshifts, it is harder for stars to keep the IGM ionized because the gas density 
is higher and there are fewer massive halos forming stars. If small halos are suppressed, the 
remaining high-mass halos have a much harder time keeping the universe ionized. It is 
interesting to note that the universe cannot be reionized at 2; = 8 — 10 if only larger halos 
(M/i > IO^Mq) are producing ionizing photons that escape into the IGM. At 2; = 9, /* must 
always be above 0.1 for all cases shown in order for the universe to be reionized. At z = 7, /* 
can be very low. Alternately, one can set equations [T7] and [12] equal and constrain directly 
the value of fescf*QLyc/CiGAi- Any values greater than this would be able to reionize the 
universe. This is shown in the bottom panel of Figure O 



Observations of high redshift galaxies at z = 7—10 (IRobertson et al.KOlOuBouwens et al 



2OIOI ) show that currently observed galaxies, along with an escape fraction that is compatible 
with lower redshift observations (/esc ~ 0.2), are not able to reionize the universe. Therefore, 
the faint end slope of the luminosity function must be very steep to create a greater contribu- 
tion to reionization from high-redshift small galaxies, the escape fraction at high redshifts was 



clumping factor of the IG 


M was low (Bunker et al. 2004 


2010: Gneiden 


2008; 


Gonzalez et al. 


lOlol: Lehnert & Bremei 


2OO3I McLure et al. 


2OI0I; 


Oesch et al. 


2010a 


b 


Oi 


chi et al. 2009: 


Richard et al. 20081 Stiavelli. Fall. & Panagia 


2004: 


Yan & Windhorst 


2004 


). In addition. 



observations of Lya absorption toward quasars and the UV luminosity function suggest that 
/esc is all owed to be much lower if lo w-mass galaxies are allowed to be sources of ionizing 
photons ( ISrbinovsky fc Wyithd 120081 ). This is consistent with our findings; small galaxies 
need to be included and their star formation not suppressed to permit /esc to be sufficiently 
high to ionize the universe. Low values of /esc, consistent with the low redshift universe, 
would require high, and in some cases unreasonably high, values of /*. Therefore, the escape 
fraction was almost certainly higher in the past than it is today. 



6. CONCLUSIONS 

We have explored how the internal properties of galaxies can affect the amount of 
escaping ionizing radiation. The properties of clumps within the galaxy have the strongest 
effect on /esc- When the covering factor is much less than one, only a few, dense clumps 
exist. These clumps only stop a small fraction of the ionizing radiation, and therefore, the 
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escape fraction will be large. As the covering factor increases to one, the escape fraction will 
fall, eventually to become smaller than the case with a non-clumped galaxy. (The exception 
to this is if the non-clumped galaxy has an escape fraction of about one. In this case the 
addition of clumps will only cause the escape fraction to fall around a covering factor of one.) 
This indicates that the escape fraction is very sensitive to the way the ISM is distributed, 
causing the escape fraction to vary from 1 to in some cases. 

The escape fraction depends on the star formation efficiency and the population of stars 
formed at the center of the galaxy. As the number of ionizing photons increases, the critical 
angle at which photons can escape the galaxy will also increase, which will have a direct 
effect on the escape fraction. Since disks were more likely more dense at higher redshifts, 
the escape fractions will be lower, unless a change in the star formation efficiency or stellar 
population would increase the number of ionizing photons. Therefore, the escape fraction 
depends directly on the number of ionizing photons and the distribution of hydrogen. These 
variations can be extreme and can help to explain why some observations see large differences 
in the escape fraction from galaxy to galaxy. Other factors, such as the mass of the galaxy 
and the redshift, affect the escape fraction indirectly (as the number of ionizing photons and 
distribution of hydrogen will change with mass or redshift). 

At high redshifts, small galaxies are probably needed to help complete reionization, 
consistent with recent observations. In addition, low values of the escape fraction, /esc ~ 0.02, 
similar to what is seen at low redshifts, are probably insufficient to allow reionization to be 
completed. Otherwise, they would require very high values of the star formation efficiency. 

We acknowledge support from the University of Colorado Astrophysical Theory Program 
through grants from NASA (NNX07AG77G) and NSF (AST07-07474). 
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Fig. 1. — Escape fraction of ionizing photons out of the disk as a function of the clump 
volume filling factor fv for various values of the clumping factor C [top panel) and log(C) 
for various values of fv {bottom panel). Shown for a IO^Mq halo at Zf = 10, with /* = 0.5 
and Population III stars with a Larson mass spectrum. The fraction of matter incorporated 
into the disk, nid, scales with the escape fraction. 
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Fig. 2. — Effect of large clumps on escape fraction: /esc out of the disk is shown for a galaxy 
in a IO^Mq halo, with zj = 10, = 0.1, and a Pop III Larson initial mass spectrum, with 
= 0.17. If the clumps are very large and fv is very low, there are cases where /esc is 
larger than the no-clump case. In the top left corner, the results are averaged over ten runs 
to reduce noise. In the remaining three panels, the distributions are shown for each run. In 
each case, the black points represent a case with no clumps. Galaxy properties are the same 
as in Figure [TJ 
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Fig. 3. — Top: The /esc out of the disk is shown for a galaxy in a IO^Mq halo, with 
Zf = 10, /* = 0.1, and a Population III Larson initial mass spectrum. The left-most 
solid vertical line represents the value of the volume filling factor fy needed for a photon 
transversing the longest path length to pass through an average of one clump for a galaxy 
with rrid = 0.17, while the right-most solid vertical line represents the value of fv needed for 
a photon transversing the shortest path length to pass through an average of one clump for 
a galaxy with rud = 0.17. Therefore, to the right of this line, the covering factor is greater 
than one, and all path lengths intersect a clump, while to the left, there are clump-free path 
lengths out of the galaxy. For a galaxy with = 0.01 the dashed vertical line represents 
the value of fv needed for a photon transversing the shortest path length to pass through an 
average of one clump for a galaxy. The value of fv needed for a photon traversing the longest 
path overlaps with the case where rrid = 0.17, because as changes, the scale height of the 
galaxy changes, but the radius remains the same. Bottom: The dependence on the escape 
fraction as rud is varied. The same population is shown, with C = 1 for the solid lines and 
C = 1000 for the dashed lines. The results are averaged over ten runs to reduce noise. 
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Fig. 4. — Distribution of /esc for each run shown in the top panel of Figure[3], with = 0.17 
in the top panel and = 0.01 in the bottom panel. 
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Fig. 5. — Escape fraction when clumps are distributed inhomogeneously throughout the 
galaxy. Two regions are defined, the region within 100 pc of star formation, which has a 
volume filling factor of /v,near, and the region outside of this, where the volume filling factor 
is /v,far- Shown for a IO^Mq halo at Zf = 10, with = 0.5 and Pop III stars with a Larson 
mass spectrum and = 0.17. Clumps are 10^^ pc across. Please see the electronic edition 
for a color version of this plot. 
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Fig. 6. — The fesc for the disk is shown for stars of varied masses and metaUicities and various 
values of fy with = 0.5 [top panel) and various values of /* for fy = 0.1 [bottom panel). 
Very high values of (0.9) approach a case in which all ionizing photons are escaping. Both 
are shown for a IO^Mq halo at Zf = 10. In each case, the highest line is for Population 
III Larson, followed by Population II Larson, Population III Salpeter, and Population II 
Salpeter. For = 0.3, /esc ~ for both Salpeter cases. Please see the electronic edition for 
a color version of this plot. 
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Fig. 7. — The /esc for the disk is shown for various values of Zf. Shown for a halo with 
Pop III stars with a Larson mass spectrum and /* = 0.5. Please see the electronic edition 
for a color version of this plot. 
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Fig. 8. — Top: The /esc needed from a population of galaxies with various values of to 
reionize the universe at 2; = 7 or 10. If the required /esc lies above 1 (dashed line), the popula- 
tion cannot reionize the universe. Bottom: The product fcscf*QLyc/CiG]\i, which is directly 
constrained by reionization. The region above each line has a value of fescf*QLyc/CiGM 
large enough to reionize the universe. Please see the electronic edition for a color version of 
this plot. 



